Crosstalk in Multiple-Beam Waveguides 

By DETLEF GLOGE 

(Manuscript received August 6, 1969) 

Crosstalk limits the number of communication channels which are 
spatially resolvable at the end of a beam waveguide. The main sources of 
crosstalk are scattering and distortion by the focusers. A careful study of 
high quality front surface mirrors led to the results of this -paper. The best 
choice seems to be a periscopic guide made of dielectric mirrors when used 
with gaussian beams in a particular mode of multiple transmission. We 
give a closed description for the expected power profile of a gaussian beam 
that has passed such a guide and an approximate formula for the mutual 
crosstalk between several such beams. 

I. INTRODUCTION 

Arranging many optical channels spatially resolved in the same 
waveguide is a simple means for high capacity transmission. 1 ' 2 All 
channels can be modulated in the same frequency band as long as the 
crosstalk is kept below a certain threshold. One source of crosstalk 
is the inevitable scattering from the focusing and directing elements. 3 
Diffraction from the ideal beams is negligible. 4 Yet we shall see that 
diffraction must be considered once the beams are distorted by the 
focusers. 

In all likelihood, these focusers would consist of mirrors because 
lenses of the size needed are apt to have imperfections within the ma- 
terial. The scattering characteristics of high quality front surface mirrors 
and lenses of the best quality have been measured recently. 3-6 A com- 
parison shows that the lens scattering was about one order of magnitude 
larger. Directional changes can easily be accomplished by using peri- 
scopic mirror arrangements of the kind shown in Fig. I. 7 Neglecting 
aberrations, we consider only the first order focusing effect of these 
periscopes, which is that of thin convex lenses. 

Two methods of multiple channel transmission have been suggested. 2 
We discuss these two basic methods with respect to their susceptibility 
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Fig. 1 — Sketch of a periscopic waveguide. 

to crosstalk when mirrors of the kind measured in Ref. 5 are used in 
the waveguide. 



II. MULTIPLE CHANNEL TRANSMISSION 

The useful cross section of the beam guide is limited by the size of 
the periscopic mirrors. Without unreasonable effort, mirrors of good 
optical quality can be made 20 to 30 cm in diameter at the most. The 
bundle of beams must clear the mirror edges by a wide margin at all 
times to guarantee safe operation. This implies that diffraction cross- 
talk caused by the mirror edges is negligible. Tolerances which would 
allow for controlled diffraction, as suggested by Ref. 4, seem unreason- 
ably tight. Thus we arrive at a radius R of about 10 cm for the useful 
cross section. The spacing D of the focusers is limited by the terrain 
and the cost of the straight sections in between. It will most likely be of 
the order of 100 m. For optimum conditions, the effective focal length 
of the focusers should be half their spacing, although some deviation 
can be tolerated in this respect. 

Consider the waveguide as a periodic lens system which images an 
array of transmitters into a similar array of detectors. This is basically 
what the imaging method (in Fig. 2a) does. Diffraction effects can be 
minimized if every transmitter radiates a coherent gaussian beam. As 
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these beams propagate in the guide, their sizes change periodically from 
the fairly small transmitter spot size to a size close to the cross section 
of the guide. 

This periodic change is avoided by the grouping method sketched in 
Fig. 2b. The beams arrange in groups and open up to the fundamental 
mode radius 



w = 



D\ 

IT 



(1) 



before they enter the guide. They keep very close to this radius through- 
out the transmission. Figure 2b shows the grouping method for two 
groups containing two beams each. Special collector lenses single out 
the groups at the end and focus the beams well separated on the de- 
tector array. For a better understanding of this detection system, 
consider the focal length / of the collector lenses to be short compared 
to the distance D between a collector lens and the preceding focuser. 
In the plane of this focuser, all groups of beams form overlapping 
patterns. Every collector lens selects the pattern of its group and images 
it into the detector plane scaled down by a factor j/D. Consequently, 
the detector array of every group is confined to a circular area with a 
radius Rf/D. 

The density of resolvable beams in the system is determined by beam 
distortion and scattering rather than the spread of the ideal beams. 
The distortion of the beam profile determines the receiver size required 




Fig. 2 — Schemes for spatially resolved transmission (a) the imaging method 
(b) the grouping method. 
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to secure safe reception. One would like to make the detectors as small 
as possible in order to minimize the scattering collected from adjacent 
channels. For further reduction of the crosstalk, one has to increase 
the spacing between the detectors. 

Let us assume that the center-to-center spacing is s for the collector 
lenses and a- for the detectors in a group. In this case we would have 
wR 2 /s 2 groups and irR 2 f/D 2 <T 2 beams per group if the guide were per- 
fectly confocal. If we allow a slight tolerance for the spacing of the 
focusers, Ref. 2 shows that the groups belonging to off-center collector 
lenses cannot be completely filled. For this reason the total number of 
channels is only half the theoretical maximum, namely 

2 s a U 

Rather than considering the detector plane, let us look at the distri- 
bution every group has at the focuser preceding a collector lens. This 
way our results become independent of the focal length / of the collector 
lenses and a function of the beam waveguide only. In the plane of the 
last focuser the beam spacing is Da/j. As the beams have the same 
width there as at the collector lenses, it seems reasonable to set 

• = Da/f. (3) 

Inserting this into equation (2) yields 

N = | (f)'. (4) 

In the following let us assume that the detectors are simple quantum 
counters, have a circular area, and have a radius df/D. Transforming 
this back to the last focuser, we find a circular area of radius d suscepti- 
ble to crosstalk around every beam. In Section II we evaluate scattering 
and distortion of a single beam in the plane of the last focuser for the 
case of a waveguide of n mirror periscopes. Since we use direct de- 
tection, we may neglect phase front distortion. The case of heterodyne 
reception is briefly discussed in the appendix. The results are veri- 
similar to those of the direct detectors. 

II. DISTORTION AND SCATTERING 

Both distortion and scattering are a consequence of irregularities 
on the mirror surfaces. The distortion originates from smooth imper- 
fections extending over areas comparable to the beam cross section 
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while scattering is caused by a surface roughness correlated over dis- 
tances much smaller than the beam diameter. Both irregularities are 
part of a statistical function 8(x, y) which describes the deviation from 
the ideal surface. Taking a meaningful average over an ensemble of 
test surfaces leads to the structure function 

A(p) = ([*(*, y) - 8(x - p cos a,y - p sin a)]\ v (5) 

where p and a belong to a polar coordinate system which has the point 
(x, y) as its origin. Writing A as a function of p only implies the assump- 
tion that 5 is stationary and isotropic, which seems justified for the 
statistical properties involved. 

A light wave of wavelength X reflected off the imperfect surface 
suffers a phase front distortion 

v(x,y) = y«(*,y). < 6 > 

We neglect reflection loss which we assume to be uniform over the 
surface. For gaussian statistics 4 

(exp i[<p(x, y) - <p{x - p cos a, y - p sin a)]) av = exp \—yT a (p)J " 

(7) 

This equality will be used to calculate the power distribution p,(r) at a 
distance D from the reflecting surface. Assume that the reflected beam is 
circular, symmetric, and would have a power profile p (r) at a distance 
D if the reflection were ideal. Then, from Ref. 8, one obtains 



gx(p) = (7o(p) e.xp 



Y- A(p)] (8) 



where g x {p) and g {p) are the Hankel transforms of p,(r) and p (r), 
respectively. This Hankel transformation is defined by 

9l (p) = 2x f X p,(r)J (2Tpr/D\)rd?- (9) 



or 



Pi(r) = ^ 5 [ gi (p)Jo(2irpr/D\) P dp (10) 

where J is the Bessel function of zero order. The quantity p, (r) has to be 
understood as an average over an ensemble of equivalent surfaces. 
For an accurate confocal spacing of the periscopes, a beam and its 
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distortion in the guide reproduces itself at every second periscope. 
These periscopes only contribute to the phase front distortion in the 
detector plane, while all odd ones deteriorate the power profile as well. 
We have n periscopes with 2n surfaces, half of them contributing to 
the profile distortion. Since the imperfections of all surfaces are un- 
corrected, we may write 

g n (p) = g (p) exp I — ^-» A(p)J (11) 

for the detector plane. 

In a guide with thousands of focusers, accurate confocal spacing 
requires tight tolerances for the focal lengths and spacings. In a practical 
guide, the focusers will be kept only nearly confocal and, in general, 
will not be at positions at which previous distortions are reproduced. 
If all positions are equally probable along the guide, equation (11) can 
be adapted in the following way 8 



g(p) = Po(p) exp 



[-^-jr'A(psinf)^]. (12) 



Notice that for £ = or w, we have A(0) =0 and no change of the 
power distribution, while for £ = t/2 the profile distortion is a maximum. 

A fairly reliable functional approximation for A in the range p = 0.01 
to 1 mm was derived from scattering measurements around a test 
beam. 6 The scattering is an effect of the mirror surface roughness 
averaged over the area covered by the test beam. This average is 
equivalent to an average over an ensemble of test surfaces. As a conse- 
quence, the variance of the scattered power is small, that is, the scat- 
tered power actually measured is very close to the average power. The 
measurements were practically the same for all test surfaces. 

This is not true for the processes involved in beam distortion. In 
this case, the 5-components participating are correlated over areas 
comparable to the test beam, no averaging is accomplished by the 
measurement, and the result can be grossly different from one surface 
to the next. It is this difference between scattering and distortion 
which makes scattering measurements feasible but distortion measure- 
ments tedious and expensive. Distortion is not sufficiently described 
by an average power profile; instead one needs to know the complete 
probability distribution of the power at every point in the beam cross 
section. In this situation some grossly simplifying assumptions are 
necessary to tackle the distortion problem with the scant experimental 
evidence available. 
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We derived the functional approximation 

A(p) = Ep for p = 0.01 • • • 1 mm 



with 



E = 2.4 X 10"°m 



(13) 



(14) 



from measurements. 5 This approximation is plotted in Fig. 3. For 
larger p we have only one reference point: the quality factor of the 
mirror, given in fractions of the green wavelength, which specifies 
5-components correlated over areas comparable to the polishing tool. 
This will be typically several centimeters. We know A decreases for 
smaller p and merges into the linear function (13) at p = 1 mm. As a 
convenient approximation, let us assume that A is linear everywhere. 
This function would correspond to about X gre e n /50 at several cm. 

Whether equation (13) is a good approximation for p < 0.01 mm 
we do not know, but this is of little interest, since components at 
these small 5 generate scattering which does not reach the next focuser, 
but is absorbed by the guide wall. 

A gaussian beam of unit power has the profile 



p (r) = — 2 exp(-2t 2 /w 2 ). 



(15) 
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Fig. 3 — The linear structure function A(p) calculated from scattering profile 
and quality tests. 
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Its Hankel transform is 

g (p) = exp (-7rVp72Z>V). (16) 

Using equations (10), (12), (13), and (16), we find the following average 
power distribution at the detector 

2 /*" ( u 2 + au\ , (2ur\ . „_> 

p (r) - w I exp v — —Mv> rfM (17) 



where 



^ J C4 W ^ /ION 

m = -=r- p and a = 64 — (18) 



The average power falling outside a circular area of radius z is 

P(z) = 2tt J p(r)r dr. (19) 



Using the identity 



we arrive at 



J,{z)z = I J (v)vdv, (20) 



P^-X-l/X-"^)^)^. (2D 

The result of the machine evaluation of equation (21) is plotted in 
Fig. 4. For a = the beam is undistorted and P(z) is gaussian. Yet 
P(z) has a tail decreasing with X/z for finite a. 

In the course of our calculation we want to know the radius z out- 
side of which a given power P can be found for a certain parameter a. 
For this purpose the function z(P, a) is plotted in Fig. 5. It can be 
approximated by the expression 

where equation (18) was inserted for a. 

The first part of equation (22) is an inverse gaussian and depicts 
the coherent beam, while the second part accounts for the incoherent 
portion. Equation (22) can be used, for example, to calculate the 
detector radius required at the end of a periscopic guide. In this case 
one will probably allow the second term in equation (22) to be about 
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Fig. 4 — The power fraction P expected outside the circular area with the nor- 
malized radius z/w for various distortion coefficients a. 



equal to the first. This would mean that the beam deterioration be- 
comes just noticeable, but not yet dominating, at the end. 

We shall find another application for equation (22) in the course of 
calculating the scattering crosstalk. In this case we require P to be so 
small that the second term in equation (22) exceeds the first even for 
moderate distortions, and equation (22) can be approximated by 



_ lQEnD 



(23) 



Notice that the only guide dimension that enters into this formula is 
the total transmission distance 



L = nD 



(24) 



from one repeater to the next. Equation (22) leads to an estimate for 
the detector sizes required and with this information and the help of 
equation (23) we can evaluate the scattering crosstalk. 
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Fig. 5 — The functional relation z(a) with Pasa parameter. 



IV. BEAM SPREAD AND CROSSTALK 

The interchannel crosstalk at the end of the multiple beam guide is a 
function of the detector size. In order to minimize the amount of light 
collected from other channels, the detectors should not be larger than 
is absolutely necessary for signal detection. A few percent of the signal 
power could even be sacrificed. The signal fraction to be detected will 
depend on the signal levels available and on the noise sources involved, 
but it is probably safe to assume that, on the average, 75 percent of the 
total signal power will be sufficient. Thus we obtain from equations (22) 
and (24) for the detector radius 



d = \w 2 In 2 + 15^16 



LE 



(25) 



In a more general sense, we may interpret d as the average radius of a 
distorted beam at the end of a guide of length L. In equation (25), w 
is the radius of the ideal gaussian beam which may vary considerably 
along the guide as, for example, in the case of the imaging method. 
For the grouping method, w is constant and given by equation (1). 

To compare both methods, let us consider a practical example of a 
waveguide with 100 m section length operating at a wavelength of 
1 A»m over a distance of 50 km. If we use the grouping method, we find 
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w = 5.65 mm from equations (1) and (25) yields a beam radius of 
8.8 mm. In the case of the imaging method, w varies about 5.65 mm 
from lens to lens, being much smaller than 5.65 mm at the detectors. 
Yet at the end of a 50 km path, the average radius of the distorted 
beams will not be smaller than 7.5 mm because of the second part of 
equation (25). This is only slightly less than the radius of the grouped 
beams. It is obvious from Fig. 2 that under these circumstances the 
imaging method loses its advantage. Actually, in this case, the imaging 
method can only accommodate the beams contained in one group of 
the grouping method. Therefore, in the following we consider only the 
grouping method. 

For the calculation of the scattering crosstalk we restrict ourselves 
to paraxial beams. Any two beams of this kind are equivalent in the 
sense that the amount of light scattered from a beam 1 into another 
beam 2 is equal to the amount scattered from 2 into 1. In the same 
way, the scattering from one beam into all others is what the beam re- 
ceives from all others. This is exactly the crosstalk we want to calcu- 
late. Thus, in order to consider the worst situation, let us select the 
center beam and calculate what it scatters into all extraneous receivers. 
To do this we have to integrate the scattered power falling into the 
detector plane, excepting the center detector and the blind area between 
all detectors. We remember that the detectors have a radius d and are 
spaced by a distance s center to center. We obtain a reasonable and 
conservative approximation if we collect all the scattering outside a 
circle with radius s/2, which is P(s/2) from equation (21), and multiply 
this by a density factor ird 2 /s 2 . Consequently, the crosstalk which the 
center beam inflicts upon, and receives from, all other beams is 

C = ^-P( s /2). (26) 

For all practical cases the tolerable C is so small that we may use the 
approximation (23) for P. Inserting equations (1), (24), and (25) into 
equation (26) we obtain 

c = 3 2^1„2 + 30.(^) 3 . (27) 

Figure 6 shows the signal to crosstalk ratio 1/C plotted in decibel 
versus the spacing s for the previous example, that is, D = 100 m and 
X = 1 nm. Also shown is the guide capacity N to be achieved by the 
grouping method in a guide of 10 cm radius. The transmission distance 
L between repeaters is the parameter. For reasons explained previously, 
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Fig. 6 — Crosstalk versus the beam density and the capacity for various re- 
peater spacings (\ = 1 /im, D = 100 m, and R = 100 mm). 



the capacity achieved by the imaging method is less, for practical sys- 
tems only about NK Figure 7 shows the useful guide radius required 
for a certain capacity at various repeater spacings if a crosstalk level 
of 23 dB is tolerable. Both Figs. 6 and 7 exhibit a system with N = 100, 
L = 50 km, R = 10 cm, D = 100 m, X = 1 am, and C = 23 dB as 
feasible but also (more or less) as a practical limit. 

The grouping method uses collector lenses in front of the detectors. 
Diffraction at these lenses must not cause excessive crosstalk even if 
the beams are badly distorted. For this reason the apertures have to be 
fairly large. If the available space is fully used, the lenses touch one 
another and are arranged as in a fly's eye lens. The lenses should be 
so large that the power at the lens edges is mainly incoherent and not 
part of the coherent, though distorted, beam. In this case, diffraction 
does not substantially increase the total power outside the signal 
beams. This requirement sets a lower limit to the beam spacing s 
which is simultaneously the diameter of the collector lenses. How far 
this limit is approached by the system depicted in Figs. 6 and 7 is a 
difficult question to answer. 

A qualitative approach is tried in Fig. 8 where the previous results 
are also summarized. Figure 8 shows the power expected outside a 
given aperture after a beam has passed a length L of periscopic wave- 
guide. The beam is supposed to start with a fundamental mode radius 
w = 5.65 mm in a guide with 100-m section length. Also shown is the 
power P(s/2) falling outside a collector lens of radius s/2 where s is 
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chosen with the help of equation (27) to guarantee a crosstalk level 
of C dB. Thus, once we have decided on the crosstalk level and the 
transmission distance, we find the radius of the collector lens and the 
power falling outside this lens from Fig. 8. For L = 50 km and C = 23 
dB, this power seems to be composed mainly of scattered light so that 
diffraction should contribute little to the overall crosstalk. 

Several discrepancies become apparent when the results of this 
paper are compared to previous publications. The power P falling out- 
side a circle with a radius z after only one reflection is obtained if we 
replace equation (12) by equation (11) and set n = 1 in the derivation 
of equation (23). In the case of a linear structure function, equations 
(11) and (12) differ by a factor ir/4, and therefore 



P = 4 .&■ 

\z 



(28) 



The same physical problem was approached on a different course in 
Ref. 5 and is expressed in equation (16) there. That result differs from 
our equation (28) by a factor of four. The reason is a factor of four 
erroneously introduced in equation (13) of that publication. 

In Ref. 3 the crosstalk of one beam into one other beam was meas- 
ured. This quantity can be calculated on the basis of tins paper. The 
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power density at a distance z from a scattering beam is 

J^dP 

2-kz dz 



Q = 



(29) 



with P{z) from equation (23). A small aperture of radius d displaced by 
a distance s from the beam center collects approximately 



c = ird 2 q(s) 



(30) 



with q from equation (29). This is the crosstalk in the second beam. 
With equations (23), (29), and (30) we obtain 



c = 



8d 2 EL 



(31) 
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The specific dimensions in the experimental arrangement of Ref. 3 
were s = 5 mm, d = 2.5 mm, and X = 0.63 /mi. The transmission 
distance was equivalent to 8.5 km of a periscopic guide. By inserting 
these data into equation (31) we obtain a crosstalk level of 19 dB. 
The measured level was 30 dB. Moreover, equation (31) suggests that 
the crosstalk decreases with the third power of the beam spacing. In 
Ref. 3, on the other hand, a decrease with the fourth power of the 
beam spacing was measured. The reason for the discrepancy is still 
under investigation. The comparison seems to indicate that the data 
used here are on the conservative side. 

Finally let us compare our results to the diffractional crosstalk which 
ideal gaussian beams experience along a wave guide. Reference 4 
considers this situation giving the following results. A guide of 10,000 
focusers, 100 m apart and 5 cm in diameter, could accommodate 16 
beams with only 60 dB crosstalk. According to our findings, scattering 
in this arrangement causes 50 dB crosstalk in one 100 m section. This 
underlines the severe limit which scattering and distortion set to multi- 
beam transmission. Improving the optical surfaces would be a worth- 
while undertaking. 

V. CONCLUSIONS 

Scattering and distortion of the beams in a beam waveguide are 
caused by surface irregularities of the focusers. There is experimental 
evidence that these irregularities can be described, as a first approxi- 
mation, by a linear structure function. Based on these findings we 
predict a beam distortion and an incoherent background radiation 
which both increase with the transmission distance L. The distortion 
makes it impossible to use a simple transmission method which images 
an array of transmitters into the detector plane. The method which 
seems more practical arranges the beams in groups and transmits them 
with ^invariable width. 

The incoherent background of scattered light around every beam 
fades off with the third power of the distance from the beam. This 
causes a crosstalk inversely proportional to the third power of the 
beam spacing. The beam density is limited by the crosstalk tolerable 
after 50 km. If we set this level at 23 dB, allow a beam bundle 20 cm in 
diameter, and operate at a wavelength of 1 /im, we could accommodate 
about 100 beams. This is based on the assumption that periscopic 
focusers are used which are made of high quality front surface mirrors. 
A critical comparison with previous publications suggests that our 
results are conservative, if not pessimistic. 
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APPENDIX 

Heterodyne Detection 

One might consider heterodyne detection as a way to reduce the 
scattering received in every channel. Local oscillator beams could be 
brought in line with the signal beams, utilizing beam splitters. The local 
oscillator beams would discriminate to a certain extent against the 
incoherent background of scattered light from other channels. To 
compare this with the quantum counters discussed in the text, let us 
assume that the scattered background light is uniform in the vicinity 
of the local oscillator beam. For this case Siegman has calculated the 
IF-photocurrent noi.se of heterodyne reception. 9 He found that the 
"integrated effective detector area" of the heterodyne detector is 
equal to X 2 . 

We now calculate the "effective detector area" for our quantum 
detectors. Every detector has an area 

A, = ^ (32) 

and is preceded by an aperture with the area 

A 2 = tts 2 . (33) 

The distance between aperture and detector is /. From Ref. 10 we find 
the "effective detector area" for this arrangement to be 

A t A a *W 



f D 2 

Inserting equation (1) we obtain 



(34) 



AiA 2 s d 2 /qi-\ 

-j- o-rX. (35) 

If we had sd = w 2 , the quantum counter would discriminate as well 
against scattering as the heterodyne detectors. 

In the case of distorted signal beams, however, both schemes cannot 
recover the full signal power. What is important in this case is the 
ratio of signal to background light collected in the respective cases. 
For this reason, the quantum counter can be equivalent to the hetero- 
dyne detector even if s and d are slightly larger than w. However, for 
reasons explained in Section IV, s will be considerably larger than w to 
avoid diffractional crosstalk. The heterodyne receiver therefore surpasses 
the quantum counter. On the other hand, the complexity of the former 
probably makes up for this advantage. 
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